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ABSTRACT 

The MAGIC telescope observed the region around the distant blazar 3C 66A for 54.2 hr in 2007 
August-December. The observations resuhed in the discovery of a 7-ray source centered at ce- 
lestial coordinates R.A. = 2'>23'°12'' and decl.= 43°0.'7 (MAGIC J0223+430), coinciding with 
the nearby radio galaxy 3C 66B. A possible association of the excess -with the blazar 3C 66A 
is discussed. The energy spectrum of MAGIC J0223-I-430 follows a power law with a normal- 
ization of (1.7±0.3stat ±0.6syst) X 10"" TeV^i cm^^ g-^ at 300 GeV and a photon index F = 

-3.10±0.31stat±0. 

Subject headings: gamma rays: observations — BL Lacertae objects: individual (3C 66A) — galaxies: 
individual (3C 66B) — ISM: individual (MAGIC J0223+430) 
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1. INTRODUCTION 

As of today, there are 23 known extragalactic very 
high energy (VHE, defined here as £^ > 100 GeV) 7-ray 
sources. All of them are active galactic nuclei (AGNs) 
with relativistic jets. With the exception of the ra- 
dio galaxy M 87 all detected sources are blazars, whose 
jets (characterized by a bulk Lorentz factor F ^ 20) 
point, within a small angle {6 ~ l/F), to the ob- 
server. The spectral energy distribution (SED, loga- 
rithm of the observed energy density versus logarithm 
of the photon energy) of AGNs shows typically a two- 
bump structure. The lower-energy bump originates 
from synchrotron radiation of relativistic electrons spi- 
raling in the magnetic field of the jet. For the origin 
of the high-frequency bump, various models have been 
proposed, the most popular invoking inverse Compton 
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scattering of ambient photons. There have been sev- 
eral suggestions for the origin of the low-frequency seed 
photons that are up-scattered to 7-ray energies: they 
may be produced within the jet by synchrotron radia- 
tion (synchrotron s elf-Comp ton or SSC m echanism, e.g. 
iMaraschi eral]fl992 : Bloom fc Marscher Iil99& ) or come 
from outs ide the jet (external Compto n or EC mecha- 
nism, e.g. lDermer fc Schlickeise"r1ll993() . Relativistic ef- 
fects boost the observed emission as the Doppler factor 
depends on the angle to the line of sight. For sources 
with a large angle between the jet and the line of sight 
(e.g., the radio galaxy M 87), these classic inverse Comp- 
ton scenarios cannot account for the VHE 7-ray emission. 
In this case, models that d epend less critically on beam- 
ing effects are needed (e.g. Neronov fc Aharonian |[2007l : 
[Tavecchio fc Ghisellini 200M). The VHE 7-ray emis- 
sion of AGNs might also be of hadronic ori gin through 
the emission from secondary electrons (e.g. Mannheiml 
[l99l iMiicke et al.i r20031 . 

3C 66A and 3C 66B are two AGNs separated by just 
6' in the sky. 3C 66B is a large Fanaroff~Riley-I-type 
(FRI) r adio galaxy , simi lar to M 87, with a redshift of 
0.0215 (|Stull et al.lll975f l. whereas 3C 66A is a blazar 
with uncertain redshift. The often referred redshift of 
0.444 (jMiller et al.lll978l ) for 3C 66A is based on a single 
measurement of one emission line only (and the authors 
were not certain on the realness of the feature) , while in 
later obse rvations n o lines in the spectra of 3C 66A were 
reported (jFinke et al. 2008 ). Based on the marginally 
resolved host galaxy (jWurtz et al.lll996l ). a photometric 
redshift of ~ 0.321 was inferred. 

3C 66A, a promising candidate for VHE 7-ray 
emission, was observed several times with satellite- 
borne and ground-based 7-ray detectors. The EGRET 
source 3EG J0222-h4 253 was associated with 3C 66A 
(jHartman et al.l Il999f ). but the association was am- 
biguous because the error box is large enough to 
cover 3C 66B and th e nearby pulsar PSR J0218+4232 
(jVerbunt et all Il996t iKuiper et all |2000D . In the 
TeV regime the Crimean Astrophysical Observatory's 
GT-48 imaging atmospheric Cerenkov telescope has 
claimed repeated det e ctions of this source ab ove 900 GeV 
(|Neshpor et al.lll998t IStepanvan et al.ll2002f ) with a flux 
as high as (3 ± 1) x lO-^cm-^s"!. HEGRA and 
WHIPPLE reported upp er limits, F(> 6 30Ge V) < 
1.42 X 10-"cm-2s-i (jAharonian et all |2000D and 
F(> 350GeV) < 0.59 x lO-^cm-^s"! (|Horan et all 
[200l . from non-simultaneous observations. The 
STACEE solar array also provided an u pper limit of 
F(> 184GeV) < 1.2 x lO'^Ocm-^s-^ ((Bramel et al.l 
I2OO5I ). In 2008 September, the Ver itas collaborat ion re- 
ported a clear detection of 3C 66A (jSwordv 1 12008') above 
100 GeV with an integral flux on the level of 10% of the 
Crab Nebula flux. Shortly after, a high state of 3C 66A 
was also reported by the Fermi Gamma-ray Space Tele- 
scope at energies above 20MeV (Tosti 2008). 

In this paper we report the discovery of VHE 7-ray 
emission located 6.'1 away from the blazar 3C 66A and 
coinciding with the radio galaxy 3C 66B in 2007. In Sec- 
tion 2, we describe the observations and the data analysis 
chain. The results of the analysis are presented in Section 
3 and discussed in Section 4. 

2. OBSERVATIONS AND DATA ANALYSIS 



3C 66 A underwent an optical outburst in 2007 August, 
as monitored by the Tuorla blazar monitoring program. 
The outburst triggered VHE 7-ray observations of the 
source with the MAGIC telescope following the Target 
of Opportunity program, which result ed in discoveries of 
new VHE 7-ray sourc es in the past (lAlbert et al.ll2006l . 
I2007at [^shhm1l2008l) . 

MAGIC has a standard trigger threshold of 60 GeV, 
an angular resolution of ^ 0.°1 and an energ y resolu tion 
above 150 GeV of - 25% (see lAlbert eFaLl (|2008afl for 
details). 

Dat a were taken in the false-source tracking (wobble) 
mode (jFomin et al.l 11993 ) pointing alternatively to two 
different sky directions, each at 24' distance from the 
3C 66A catalog position. The zenith distance distribu- 
tion of the data extends from 13° to 35°. Observations 
were made in 2007 August, September, and December 
and lasted 54.2 hr, out of which 45.3 hr passed the qual- 
ity cuts based on the event rate after image cleaning. An 
additional cut removed the events with total charge less 
than 150 photoelectrons (phe) in order to assure a better 
background rejection. 

Just before the start of the observation campaign ~ 5% 
of the mirrors on the telescope were replaced, worsen- 
ing the optical point-spread function (PSF). As a conse- 
quence, a new calibration of the mirror alignment system 
became necessary, which took place within the observa- 
tion campaign and improved the PSF again. The sigma 
of the Gaussian PSF (40% light containment) was mea- 
sured to be 3.'0 in 2007 August 12-14, 2.'6 in 2007 August 
15-26 and 2.'1 in 2007 September and December. To take 
this into account, data were analyzed separately for each 
period and the results were combined at the end of the 
analysis chain. However, the realignment resulted in a 
mispointing, which was taken care of by a new pointing 
model (Brctz ct al. 2009) applied offline using starguider 
information (jRiegel et al.ll20C)5[ l. Considering the addi- 
tional uncertainty caused by the offline corrections, we 
estimate the systematic uncertainty of the pointing accu- 
racy to be 2' on average. Note that in the case of an opti- 
mal pointing model the systematic uncertainty is below 
2\ bei ng 1' on average (jBretz et al.l l2009l : lAlbert et all 
I2008a[) . 

The data analysis consists of se veral steps. Initially , 
a standard calibration of the data (| Albert et al.ll2008dl ) 
is performed. In the next step, an image cleaning pro- 
cedure is applied using the amplitude and timing infor- 
mation of the calibrated signals. In particular, the ar- 
rival times of the photons in core pixels (> 6 phe) are 
required to be within a time window of 4.5 ns and for 
boundary pixels (> 3 phe) within a time window of 1.5 ns 
from a neighboring core pixel. For the surviving pix- 
els of each event image parameters are calculated ( Hillai 
Il985l ). Using the good time resolution of the recorded 
signals (~ 400 ps), unique to MAGIC, the time gradient 
along the main shower axis and the time sp read of the 
shower pixels are computed (Aliu_et_ al. 2009) . Hadronic 
background suppress ion is achieved usi ng the Random 
Forest (RF) method () Albert et al.ll2008cD . where for each 
event the so-called Hadronness parameter is computed, 
based on the image and the time parameters. Moreover, 
the RF method is used for the energy estimation trained 
on a Monte Carlo simulated 7-ray sample with the same 
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Fig. 1. — Significance map for 7-like events above 150 GeV in 
the observed sky region. The green cross corresponds to the fitted 
maximum excess position of MAGIC J0223+403. The probability 
of the true source to be inside the green circles is 68.2%, 95.4%, 
and 99.7% for the inner, middle, and outer contour, respectively. 
The catalog positions of 3C 66A and 3C 66B are indicated by a 
white square and a black dot, respectively. 

zenith angle distribution as the data sample. 

3. RESULTS 

Figure [T] shows a significance map produced from the 
signal and background maps, both smoothed with a 
Gaussian of ti = 6' (corresponding to the 7-PSF), for 
photon energies between 150 GeV and 1 TeV. For the 
background rejection a loose cut in the Hadronness 
parameter is applied to keep a large number of gamma- 
like events. The center of gravity of the 7-ray emission 
is derived from Figure [1] by fitting a bell-shaped function 
of the form 



F{x, y) = A ■ exp 



{x ~ x)^ + {y - y)^ 



(1) 



for which the distribution of the excess events is assumed 
to be rotationally symmetric, i.e., (J^ = Cy = cr. The 
fit yields reconstructed coordinates of the excess center 
of R.A. = 2i^23'"12'' and decl.= 43°0.'7. The detected 
excess, which we name MAGIC J0223+430, is 6.'1 away 
from the catalog position of 3C 66A, while the distance 
to 3C 66B is l.'l. 

In order to estimate the statistical uncertainty of the 
reconstructed position, we simulated 10** sky maps with 
the same number of background and excess events as in 
the data. The excess position in the sky maps was fitted 
and the distance to the simulated source position calcu- 
lated. From the histogram of the distances we obtained 
probabilities for an offset between the true source and 
the fit to the excess. The probabilities shown in Figure [T] 
by the green contours correspond to 68.2%, 95.4%, and 
99.7% for the inner, middle, and outer contour, respec- 
tively. Using this study we find that the measured excess 
coincides with the catalog position of 3C 66B. The origin 
of the emission from 3C 66A can be statistically excluded 
with a probability of 95.6%. Adding linearly the system- 
atic uncertainty of the pointing of the data set (2', see 
above), i.e., shifting the excess position by 2' toward the 
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Fig. 2. — |Alpha| distribution after all cuts evaluated with re- 
spect to the position of MAGIC J0223-I-430. A 7-ray excess with 
a significance of 6.0 cr is found, which corresponds to a post-trial 
significance of 5.4 cr. 



catalog position of 3C 66A, the exclusion probability is 
85.4%. 

To calculate the significance of the detection, an 
|Alpha| distribution was produced, where Alpha is the 
angle between the major axis of the shower image el- 
lipse and the source position in the camera. For the cal- 
culation of the source-dependent image parameters we 
considered the fitted position of the excess. Background 
rejection was achieved by a cut in Hadronness, which 
was optimized using Crab Nebula data taken in similar 
conditions and diluted to 5% of its real flux. The cut 
in I Alpha I that defines the signal region was also opti- 
mized in the same way. The | Alpha | and Hadronness 
cuts together have an efficiency of 40% in keeping Monte 
Carlo simulated 7 events, and result in an energy thresh- 
old of approximately 230 GeV. ^ A signal of 6.0 cr signifi- 
cance (pre-trial) was found (see Figure [2]). We estimated 
the number of trials of the signal search by projecting 
the 7-ray acceptance of the camera into the field of view 
of the observations, and defined the search region where 
the 7-ray acceptance after cuts is larger than 50%. In 
this way, we obtained an area of 2.18 deg^. Given that 
the 68% containment radius for 7-rays from a point-like 
source is 0.°152, we calculated the number of indepen- 
dent trials to be 30. 

Figure [3] shows the light curve of MAGIC J0223-H430 
together with the flux of 3C 66 A in optical wavelengths. 
As we integrate over 7-ray events from a wide sky re- 
gion (~ 0.07 deg^), we cannot exclude that 3C 66A con- 
tributes to the measured signal. The integral flux above 
150 GeV corresponds to (7.3 ± 1.5) x 10"^^ ^m-a g-i 
(2.2% of the Crab Nebula flux) and is the lowest ever 
detected by MAGIC. The 7-ray light curve is consistent 
with a constant flux within statistical errors. These er- 
rors, however, are large, and some variability of the signal 
cannot be excluded. 

For the energy spectrum of MAGIC J0223-f 430, loose 
cuts are made to keep the 7-ray acceptance high. 
The differential energy spectrum was unfolded usin g 
the Tikhonov unfolding te chnique (jTikhonov Sz Arsenini 
[T979I: [Albert et al.]|2007bD and is shown in Fig. H The 

^ Defined as the peak of the distribution of Monte Carlo gener- 
ated gamma-ray events after all cuts. 



4 



20 



















J, 






MAGIC 


J0223t430 






I 





















































1'' 

I 10 

i 8 

to 

SI 

d: 6 















I 




I 


3C 66A (KV; 






I J I 






I 


I 

■ I 










■ I . 


in* 


'l . I 


I 


1 

I 















54320 54340 54360 54380 54400 54420 54440 

time [MJD] 

Fig. 3.— Light curve of MAGIC J0223+430. Upper panel: 
MAGIC integral flux above 150 GeV in bins of 3 days (except for 
periods where the sampling was coarser). The gray dashed line 
indicates the average 7-ray flux. Lower panel: optical light curve 
of 3C 66A as measured by the KVA telescope. During the MAGIC 
observations 3C 66A was very bright at optical wavelengths vary- 
ing from 6 mjy to 12 mjy in the ij-band (the baseline flux in the 
historical data being ~ 6 mjy). In the same period the optical flux 
of 3C 66B remained constant, which is a typical behavior for large 
radio galaxies. 



spectrum can be well fitted by a power law which gives 
a differential flux (TeV^^ cm^^ s^^) of: 



dN 



dEdAdt 



= (1.7±0.3)xlO""(i;/300GeV) 



-3.1±0.3 



(2) 



The quoted errors are statistical only. The systematic 
uncertainty is estimated to be 35% in the flux level an d 
0.2 in the power law photon index (| Albert et al.ll2008al ). 
As in the case of the light curve, we cannot exclude that 
3C 66A contributes to the measured signal. Thus, the 
spectrum shown in Figure |4] represents a combined 7-ray 
spectrum from the observed region. 

4. DISCUSSION AND CONCLUSIONS 

A new VHE 7-ray source MAGIC J0223-h430 was de- 
tected in 2007 August to December. Given the position 
of the excess measured by MAGIC above 150 GeV, the 
source of the 7-rays is most likely 3C 66B. The VHE 
7-ray flux was found to be on the level of 2.2% Crab 
Nebula flux and was constant during the observations. 
The differential spectrum of MAGIC J0223-f 430 has a 
photon spectral index of F = 3.10 ± 0.31 and extends up 
to '--^ 2 TeV. I n vie w of the recent detection of 3C 66A at 
VHE 7-rays (S wordy •2008), we note that if 3C 66 A was 
emitting 7-rays in 2007 August to December then its flux 
was at a significantly lower level than in 2008. We also 
note that we c annot exclude the scenar i o sugg ested in a 
recent work by iTavecchio fc Ghisellini I (|2008b( ) that the 
observed spectrum would be a combination of emission 
from 3C 66B (dominating at energies above 150 GeV) 
and blazar 3C 66 A (at lower energies). 

In the unlikely case, excluded with probability 85.4%, 
that the total signal and observed spectrum presented in 
this paper originates from 3C 66A, the redshift of the 
source is likely to be significantly lower than previously 
assumed. Due to the energy-dependent absorption of 
VHE 7-rays with low-energy photons of the extragalactic 
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Fig. 4.— Differential energy spectrum of MAGIC J0223-I-430. 
The flt to the data is shown by the solid gray line and the fit 
parameters are listed in the inset. No correction for the 7 — 7 
attenuatio n due to the EBL has been made. The Crab Nebula 
spectrum IIAlbert et al.l2008al ) is also shown as a reference (dashed 
gray line). 



background fEBL. ' Gould fc Schreder |[T96l . the VHE 7- 
ray flux of distant sources is significantly suppressed. We 
investigated the measu red spectrum by MAG IC follow- 
ing the prescription of iMazin &: Rauel (|2007f ). and de- 
rived a redshift upper limit of the source to be z < 0.17 
(z < 0.24) under the assumption that the intrinsic energy 
spectrum cannot be harder than F = 1.5 ( F = 0.666). 
This assumption of F > 1.5 is based on p article acceler- 
ation arguments (|Aharo nian et al.l l2006af ). and the fact 
that none of the sources in the EGRET energy band (not 
affected by the EBL) have shown a harder spectrum. 
The latter assumption of F > 0.666 can be considered 
as an extreme case of the spectrum hardness, suggesting 
a monochromatic spectrum of electrons whe n interact- 
ing w ith a soft photon target field (Katarz vnski et all 
l2C)06ft .^ If z > 0.24 for 3C 66A, an affernative expla- 
nation for a hard intrinsic spectrum at energies above 
100 GeV can be given if 7-rays are passing through a 
narrow band of optical-infrared photons in the vicinity 
of the blazar. Such narrow radiation fields can produce 
arbitrarily hard intrins ic spectra b y absorbing specific 
energies of 7-rays (Aharonian et al.]|2008) . We also note 
that, in this case, the intrinsic VHE luminosity of 3C 66A 
should exceed 10^^ ergs" \_ which is an unusually large 
value for a BL Lac object ([Waancr 2008) , also in view of 
its spectral characteristics ( Persic &: De Angelis 1 12008*). 

3C 66B is a FRI radio galaxy simila r to M 87, which 
has b e en detected to emit VHE 7-rays () Aharonian et al.l 
[200I lAharonian eTall l2006bD . Since the distance 
of 3C 66B is 85.5 Mpc, its intrinsic VHE luminosity 
would be two to eight times higher than the one of 
M 87 (22.5 Mpc) given the reported variabil ity of M 87 
(|Aharonian et al.ll2006bHAlbert et al.ll2008b[) . 

As in the case of M 87, there would be several pos- 
sibilities for the region responsible of the TeV radia- 
tion i n 3C 66B: the vicinity of th e supermassive black 
hole (Neronov fc Aharonian I [20071 ). the unresolved base 
of the jet (in analogy with blazar emission models; 

^ See also lStecker et all I I2007I) for more detailed calculations. 
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iTavecchio fc Ghisellini |[2008a[ ) and the resolved jet. Un- 
like for M 87, we do not observe significant variability in 
the VHE 7-ray flux and therefore we have no constrains 
on the size of the emission region. However, as the angle 
to line of sight is eve n larger than in M 87 (M 87: 19° 
iPerlman et all (j2003f ): 3C 66B: 45°, iGiovannini et all 
(|2001h ) the resolved jet seems an unlikely site of the emis- 
sion. On the other hand, the unresolved base of the jet 
seems a likely candidate for the emission site as it could 
point with a smaller angle to the line of sight. If the 
viewing angle was small, blazar-like emission mechanisms 
cannot be excluded. The orbital motion of 3C 66B shows 
evidence for a supermassive black ho le binary (SMBHB ) 
with a period of 1.05 ± 0.03 years (jSudou et al.ll2003D . 
The SMBHB would likely cause the jet to be helical, and 
the pointing direction of the unresolved jet could differ 
significantly from the direction of the resolved jet. 

Given the likely association of MAGIC J0223-H430 with 
3C 66B, our detection would establish radio galaxies as 



a new class of VHE 7- ray emitting sources. According to 
iGhisellini etlHI (|2005| ). there are eight FR I radio galax- 
ies in the SCR catalog that should have a higher 7-ray 
flux at 100 MeV than 3C 66B, but possibly many of these 
sources are rather weak in the VHE 7-ray band. Further 
observations of radio galaxies with the Fermi Gamma-ray 
Space Telescope as well as by ground-based telescopes are 
needed to further study the 7-ray emission properties of 
radio galaxies. 
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